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Abstract 

A linear relationship between the Hubble expansion parameter and the time derivative 
of the scalar field is explored in order to derive exact cosmological, attractor-like solutions, 
both in Einstein's theory and in Brans-Dicke gravity with two fiuids: a background fiuid of 
ordinary matter, together with a self-interacting scalar field accounting for the dark energy in 
the universe. A priori assumptions about the functional form of the self-interaction potential 
or about the scale factor behavior are not necessary. These are obtained as outputs of the 
assumed relationship between the Hubble parameter and the time derivative of the scalar 
field. A parametric class of scaling quintessence models given by a self-interaction potential 
of a peculiar form: a combination of exponentials with dependence on the barotropic index of 
the background fiuid, arises. Both normal quintessence described by a self-interacting scalar 
field minimally coupled to gravity and Brans-Dicke quintessence given by a non-minimally 
coupled scalar field are then analyzed and the relevance of these models for the description 
of the cosmic evolution are discussed in some detail. The stability of these solutions is also 
briefiy commented. 



1 Introduction 



Dark energy or missing energy is one of the contemporary issues the physics community 
is more interested in due, mainly, to a relatively recent (revolutionary) discovery that our 
present universe is in a stage of accelerated expansion[l]. This missing component of the ma- 
terial content of the universe is the responsible for the current stage of accelerated expansion 
and accounts for 2/3 of the total energy content of the universe, determining its destiny[2]. 
This is a new form of energy with repulsive gravity and possible implications for quantum 
theory and supersymmetry breaking. Among others, the new cosmology is characterized 
by the following features[2]: 1) Flat, critical density accelerating universe, 2) early period 
of inflation, 3) density inhomogcncitics produced from quantum fluctuations during infla- 
tion, 4) composition of 2/3 dark energy, 1/3 dark matter and 1/200 bright stars, 5) matter 
content: (29 ± 4)% cold dark matter, (4 ± 1)% baryons and ~ 0.3% neutrinos. Besides, 
there is evidence (recent observation of SN1997ff at redshift z = 1.7) that the present stage 
of accelerated expansion was preceded by an early period of decelerated expansion[2]. Any 
theoretical framework intended to describe our present universe should be compatible with 
most of these requirements.^ 

A self-interacting, slowly varying scalar field, most often called quintessence, has been 
meant to account for the dark energy component. In the simplest models of quintessence, 
minimally coupled scalar fields are considered. However, extended quintessence models have 
also been considered[3, 4, 5]. In these models the scalar field that accounts for the dark 
energy component is non-minimally coupled to gravity and, besides, is a part of the gravity 
sector. Brans-Dicke theory is a prototype for this kind. Motivated by the possibility to have 
superaccelerated expansion and, correspondingly, superquintessence solutions, this type of 
theories is now being considered (see, for instance, Ref . [5] and references therein) . 

A variety of self-interaction potentials for the quintessence field has been studied. Among 
them, the simplest exponential potential (a single exponential) model is unacceptable because 
it can not produce the transition from subdominant to dominant energy density[6]. Other 
combinations of exponential potentials have been also studied[7, 8, 9, 10, 11]. Combinations 
of exponentials are interesting alternatives since these arise in more fundamental (particle) 
contexts: supergravity and superstring[12], where these types of potentials appear after 
dimensional reduction. 

In most cases the occurrence of a self-interaction potential for the scalar field makes diffi- 
cult to solve analytically the field equations, although some techniques for deriving solutions 
have been developed. In Ref. [13], for instance, the form of the scale factor is given a priori 
and, consequently, the self-interaction potential can be found that obeys the field equations. 
This method has been repeatedly used[14, 15] as well as interesting applications of it[16]. 
However there are cases when exact solutions can be found once the form of the potential is 
given[10, 11]. In other cases some suitable relationship between the self-interaction poten- 
tial and the scalar field kinetic energy is assumed [17]. The assumption that the scalar field 

^Most part of these features can be considered as observational evidence. Early inflation, on the contrary, 
is just a theoretical model that is consistent with observations. 
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energy density scales as an exact power of the scale factor has also been used for deriving 
solutions [18]. 

In case we deal with a non-minimally coupled scalar field that accounts for the dark 
energy, the situation with finding exact solutions is even more complicated. It seems that 
numerical integration is the "simplest" way to treat this problem. In Ref. [5], for instance, 
numerical solutions are explored in extended quintessence models. 

In this paper we explore a linear relationship between the Hubble expansion parameter 
and the time derivative of the scalar field to derive exact cosmological scaling solutions^ both 
in Einstein's gravity and in Brans-Dicke theory with two fiuids: a barotropic perfect fluid 
of ordinary matter, together with a self-interacting scalar field fluid accounting for the dark 
energy in the universe. 

A similar (but more general) relationship (f) oc H + aL{(f)) has been used in Ref. [20] to 
find general solutions in + 1-dimensional theories of gravity with a self- interacting scalar 
field and also to derive general solutions when one deals with two scalars: a self-interacting 
scalar field with exponential potential plus a free scalar field. ^ 

We show that the assumed linear relationship between the Hubble parameter and the 
time derivative of the scalar field is suggested by an implicit symmetry of the field equations 
and it can lead always to attractor scaling solutions. In this regard, we point out that it is not 
necessary to make any a priori assumptions about the functional form of the self-interaction 
potential or about the scale factor behavior. These are obtained as outputs of the assumed 
relationship once one integrates the field equations explicitly. As a result a peculiar class 
of potentials arise: a combination of exponentials with some dependence on the barotropic 
index of the ordinary matter fiuid. This kind of potential is of interest for models of dark 
energy since, as it has been pointed out in Ref. [12] , a single exponential potential can lead to 
one of the following scahng attractor solutions: 1) a late time attractor where the scalar field 
mimics the evolution of the background barotropic fluid, and 2) the late time attractor is the 
scalar fleld dominated solution. In consequence, a combination of exponentials shoukl allow 
for a scenario in which the universe evolves through a radiation-matter regime (attractor 1) 
and at some recent epoch evolves into the scalar fleld dominated regime (attractor 2). 

In this paper we will be concerned with flat Priedmann-Robertson- Walker (FRW) cos- 
mologies with the line element given by: 

ds^ = -dt^ + a{tfSikdx'dx\ (1) 

where the indexes i. A; = 1, 2, 3 and a{t) is the scale factor. We use the system of units in 
which 8ttG = c = 1. 

The paper has been organized in the following way. In section 2 we explain the details 
of the method by applying it to derive solutions in Einstein's theory with two fluids: a 
background fluid of ordinary matter and a fluid of a self-interacting scalar fleld that accounts 

^By scaling we mean cosmological models where the dimensionless density parameter for each component 
in the energy content is a constant: Cli = pi/pcrii = constant [19]. 

method where a linear relationship between the field variables and/or their derivatives is assumed has 
been also used in [21] to derive 4d Poincare invariant solutions in thick brane contexts. 
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for the dark energy component in the universe. We study also, the resulting parametric class 
of scaling solutions (including the space of parameters) and we briefly comment on their 
stability. In particular, we make clear the point that the assumed linear relationship between 
the Hubble expansion parameter and the time derivative of the scalar field, acts just like 
a selection principle that preserves only the (relevant) attractor-like solutions. Section 3 is 
devoted to the derivation of a parametric class of scahng solutions in Brans-Dicke gravity 
with a background fluid and a self-interacting Brans-Dicke (BD) field, by assuming the same 
linear relationship between the Hubble expansion parameter and the time derivative of the 
BD field. We apply the conformal technique to simplify the mathematical handling of the 
problem. The stabihty of the solutions found is also commented in this section, while their 
relevance for the description of the cosmic evolution will be discussed in some detail in section 
4. Some conclusions are given in the final section 5. 



2 The method 

In this section we explain the details of the method we use for deriving solutions in gravity 
coupled to self-interacting scalars. We will derive solutions in Einstein's gravity with two 
fiuids: a background barotropic fiuid of ordinary matter and a self-interacting scalar field 
fiuid. The field equations are the following: 

3H^ ^Pm + + V, (2) 
2H + = (1 - 7)p^ - + V, (3) 

4> + 3H4>^ -V, (4) 

where V = ^(0) is the self-interaction potential, 7 is the barotropic index of the background 
fluid and H = - is the Hubble expansion parameter. The dot accounts for derivative in 
respect to the comoving time t meanwhile the comma denotes derivative in respect to the 
scalar field 0. The energy density of the ordinary matter (cold dark matter plus baryons 
and/or radiation) is related with the scale factor through — po^^a~^'^, where po,7 is an 
integration constant coming from integrating the conservation equation. Let us combine 
equations (2) and (3) to obtain 

H + 3H^^ '^Pm + V. (5) 

An implicit symmetry of the left hand side (LHS) of equations (4) and (5) is evident under 
the change H k(f). 1. e., if one assumes a linear relationship between the Hubble expansion 
parameter and the time derivative of the scalar field; 

H^k^,^a^ e^'t', (6) 



3 



where A; is a constant parameter, the LHS of equations (4) and (5) coincide up to the factor 
k. We obtain a differential equation for determining the functional form of the potential V: 



Explicit integration yields the following potential which is a combination of exponentials: 

V = ^oe-^/^ + ^|^Po,.e-^^-^ (8) 

where .^o is an integration constant. This potential is a particular case of the potential V{(f)) — 
M\e'^'t> + e^-^) in Ref.[12], when a = l/k, (3 = 3A;7, and = q^^Po^ = M^- It can be 

given in terms of the scale factor if one considers equation (6): V = ^Qa~^^^ + gfc2~l2 Po.7'^~''^^- 
An interesting feature of this potential is that it depends on the type of ordinary fluid 
which fills the universe. Otherwise, it depends on the barotropic index 7 of the matter fluid. 
If one introduces the following time variable dt — a)-!'^^ dr and by substituting (8) back into 
Eq. (2), one can rewrite this last equation in the following form: 

where now the dot means derivative in respect to the new time variable r. This equation 
can be integrated in quadratures: 

^37/2-1/2^2-1^^ 

= T + To, (10) 



where A — ^^^_2 Po,ii B = ^z^j and Tq is an integration constant. 
2.1 The class of solutions 

The integral in the LHS of Eq. (10) can be explicitly taken to yield the following expression 
for the scale factor as function of the time variable r; 

a(T) = ao{sinh[/i(T + To)]y''/^''''^~'\ (11) 

where Uq = (A/5)'='/(3'='7-i) and // = {3k^-f - l)y/B/k'^. Then, according to Eq. (6) the 
scalar field can be easily computed to yield, 

1 2k 

(/)(t) = -lnao + ^p;^lnsinh[/x(T + To)]. (12) 

The Hubble expansion parameter can be computed by the knowledge of the scale factor 
(11): 



^^^^ ^ J^|^«r)-^/^'=%oth[/.(r + ro)]. (13) 



For the purpose of observational testing of the solutions it is useful to look for another 
magnitudes of astrophysical interest. Among them, the matter density parameter: 



n„ = pj3H' = ^^2__|cosh[/.(T + To)]}-^ (14) 
the scalar field energy density and scalar field density parameter: 

%={t^ = (15) 

respectively. The second equation in (15) is just another writing for the field equation (2). 
The scalar field state parameter can be written in the following form: 

\4^^-v 1 .... 

meanwhile the deceleration parameter q — —{1 + H/H'^): 

g = -1 + ^ + f o^. (17) 

While deriving equations (15)-(17) we have used the field equations (2), (3), (4) and 
their combinations. Assuming the linear relationship (6) between the Hubble parameter and 
the time derivative of the scalar field means that we are introducing a new free parameter, 
however, as it will be shown latter on in subsection 2.B, the space of parameters can be 
reduced by the combined action of choosing an appropriated normalization and considering 
some well established observational facts. For purpose of observational testing of the solution 
it is desirable to have the physical parameters of observational interest written as functions 
of redshift z. This is achieved simply by replacing in the former equations 

^^(r)^n„U) = (>l^)^^S}±Bt_, (18) 

where 2n = 87 — l//c^. While deriving this expression one has considered a{z) — a(0) /{z + 1), 
and the normalization a(0) = 1. In terms of the redshift the Hubble expansion parameter 
can be written as: 



H{z) = VaJ {z + If-r + + l)Vfe^ (19) 



meanwhile, the background energy density pm = Po,'y{z + 1)^'''. 
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2.2 The space of parameters 



We now investigate the space of parameters (7,^0,^07, k) parameterizing the solutions (the 
parameter 7 can be fixed to be 1 at present and 4/3 in the past). We consider k"^ > since 
otherwise the scalar field would be a complex one. We should differentiate two ranges. For 
< A;^ < 1/87 the scale factor evolves periodically: 



or 



~ao{sm[uiT + To)]y''/^'''^-'\ (20) 

where gq = [—A/BY^^^^^'^^~^\ However, in this case, for < 7 < 2, 3A;^7 — 1 < so one 
should have ^0 < 0. This means that the self-interaction potential itself is negative, i. e., 
the scalar field energy density is negative. Besides q > Q always and this case is not of 
observational interest. Therefore wc drop this interval in parameter k. 

For k"^ > 1/37 the equations (11)-(13) are valid in that form. To constrain further the 
space of parameters we use some observational facts. However, this point is left for further 
discussion in section 4. We now determine the constant A/B. For this purpose one uses two 
observational facts: 1) Since, at present Jlm(O) = 1/3, then one obtains from Eq. (18) that: 
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(21) 



and 2) Since radiation-matter dominated the early stages of the evolution, then one may 
require that Qrn{oo) = 1 — e, where the new parameter e << 1. This small parameter is 
necessary since a matter-dominated {flm = 1) is not stable [22]. In consequence, from the 
same Eq.(18), one obtains that: 



Po,7 = 3A(1 - e). 



Combining the above equahties one gets that: ^ — 2 — 3e. In consequence: 



n^niz) = (1 - e) 



A 

[Z + 1)2" 



(^ + l)2"-f2-3£' 



and. 



\ 



3fl-£) 



(22) 



(23) 



(24) 



One can reduce the space of parameters (7,6,^0,7;^) by requiring the following normal- 
ization: i^(0) = 1, implying that Pm(0) = Po,7 — 1- We are then left with a three-parameter 
class of solutions spanned by (7, £, k). In section 4 we will further constrain the parameters 
e and k by the observations. 



2.3 Stability of the solutions 

The stability analysis of the model given by Eqs.(2.1)-(2.3) has been performed in Ref.[22] for 
the class of exponential potentials of the kind: V — Vqc'^'^ {k — SttG = 1). There are up to 
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five fixed points (critical points). Two of them are late-time attractor solutions. One of these 
is the scalar field dominated solution: Q^j, = 1, with an afi^ective barotropic index 7^ = A^/S. 
This solution exist for < 6, and is a lite-time attractor in the presence of a background 
fluid when A^ < 87 [22]. The other attractor solution occurs for A^ > 87 and neither the 
scalar-field nor the background fluid entirely dominates the evolution.lt is a scaling solution 
where il<^ = |i . At this point we want to discuss on the implications of the linear relationship 
(2.5) assumed in order to find exact solutions. The procedure for the stability analysis [22] 
implies writing the evolution equations (2.1)-(2.8) as a plane-autonomous system. For this 
purpose, one should define new variables: 



_ ^ _ w 

^ " V6H' ^ " VSH' ^ ' 

Considering (6), the first equation implies x — Meanwhile the critical points leading 
to the former attractor solutions are xi — Xi/V^ for the scalar dominated solution and 
X2 = y^7/A2 for the scaling solution [22]. On the other hand, our potential (Eq.(8)) 
contains a combination of exponentials with exponents Ai = 1/k and A2 = 8^7. We see 
then, that assuming the linear relationship (6) means dropping of the critical points but the 
ones that can lead to attractor solutions. Otherwise, relationship (6) could be written in 
the following equivalent way: 70^^ = Then one is faced with two possibilities: either 
Jl-y = 1 ^ 7^ = gp- (scalar field dominated solution), or 7^ = 7, ^ct> — — const 
(scaling solution). In this sense one can say that the assumed relationship between H and 
(Eq. (6)) acts just like a selection principle that, among the possible critical points, preserves 
only the (relevant) attractor-like solutions. 



3 Brans-Dicke quintessence 

Now we use the former method to derive exact solutions in Brans-Dicke gravity with a 
self-interacting scalar field. As mentioned in the introduction to this paper, this is one 
in a kind of models of extended quintessence studied before in [8, 4, 5]. In Ref. [5], in 
particular, numerical solutions in extended models of quintessence (including Brans-Dicke) 
were explored. 

The Brans-Dicke field equations, in the presence of a self-interacting scalar field, are: 

3//2 + 3//(^ - = e-'^ipm + V), (26) 

2ij + 3i/2 + V3 + 02 + 2Hif + ^ _ + (27) 

+ 02 + 3//<^ = - 87)p„. + 4y - 2V'}, (28) 
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where we have introduced a new field variable </? instead of the original (j) — e'^. In these 
equations uj is the Brans-Dicke coupling parameter. 7 is the barotropic index of the fiuid 
of ordinary matter, V is the self-interaction potential and H = ^ is the Hubble expansion 
parameter. The dot accounts for derivative in respect to the comoving time t meanwhile the 
comma denotes derivative in respect to the scalar field (f. As before, the energy density of the 
ordinary matter (cold dark matter plus baryons and/or radiation) is related with the scale 
factor through = po^^a~^'*, where po,7 is an integration constant coming from integrating 
the conservation equation. 

For simplicity of the mathematical handling of the problem it is convenient to perform 
a conformal transformation of the metric: (jab = ^'^^^Qah-, where ^ = a; + 3/2 and, for 
convenience, we have rescaled the BD field (p — ^/v^. The transformed fine element is 
now: ds^ — —dP + a^Sijdx^dx^ , where the conformal comoving time and scale factor are 

related with the original (Jordan frame) ones through: dt = e and a — e~'^/^'^"^^a 

respectively. Then we will be faced with the Einstein frame (EF) formulation of BD theory. 
The relevant field equations are now the following: 

3^' ^Pm + Pi,, (29) 
2H + 3H^ = (1 - 7)p^ - P^, (30) 

^ + 3^^ = ^^p^-y', (31) 



where the EF scalar field energy density = ^tp + V and the corresponding pressure 
\ ~ 2^ 



Pip — ^"0^ ~ V- Besides, the transformed "conservation" equation is: 



^^ + {3-fH + ^'^Wi;)pm^0, (32) 

where the "exchange" function W = ^f^^ has been introduced as in Ref.[19]. After 
integration: = po^ja~^'^e~^^^ . Summing up (29) and (30) yields: 

S + SH^^'^p^ + V. (33) 

As before an implicit symmetry of the LHS of equations (31) and (33) is made explicit 
by the change ip XH. In this case Eq.(31)^ A times Eq.(33). Therefore, if we assume the 
linear relationship 

^ = A^, e^ = a^, (34) 

then the following differential equation for the self-interaction potential V is obtained: V' + 
\V — [2J^W — (2 — 7)A]pr„/2. Straightforward integration of this equation yields: 
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V{ij)^Voe-^'l' + Woe-^'^, (35) 

where S — 3^/X + ^2/3W, Vq is an integration constant and Wq — [^J^W — ^^\]po,'y/{X — 5). 

In terms of the scale factor Eq. (35) reads: V{ip{a)) = Vb«~^^ + Woa~^^. 

This potential belongs to the same class that (8) and shows the interesting feature that 
it depends on the type of ordinary fluid which fills the universe (it depends on the barotropic 
index 7 of the matter fluid). By substituting Eq. (35) back into the Friedmann equation 
(29) one obtains 



^ Aa-^^ + Ba-^\ (36) 

where ^ = (A^ - G)-f/{\'^ - ^l2/2,W\ - 37), B = 2Vo/{Q - A^), and we have considered, 

besides, that = pQ^^cT'^'^'^^^ = pQ^^a^^^. 

To solve (36) we make another change of time variable: dr — a'^^^^dt to get Eq. (36) 
integrated in quadratures: 

where ro is another integration constant. 

3.1 The class of solutions 

Straightforward integration of Eq. (37) leads to: 

^^^y{5-x)/2 ^ JA/Bsm\i[ii{r + Tq)], (38) 



with // = B/^X{5 — A)/2. For the rescaled BD scalar field we obtain: 



= {^A/Bsinh[p(r + r^)]Y'^^'^\ (39) 
In consequence the dimensionless background density parameter: 



and. 



Vt^ = l- Clrn- (41) 

As before it is useful to have the magnitudes of observational interest written in terms 
of the redshift z. In the Jordan frame we have that a{z) = a{0)/{z + 1). We choose the 

normalization a(0) = 1 ^ a(0) = 1. Then, since a = e~'^l^^^^^a — a"' — >• a = {z + 1)~^^", 
where we have introduced the following definition for the constant n: 
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After these considerations the dimensionless density parameter and the Hubble expansion 
parameter can be written in the following way: 

n (z) = ^ ^) C43A 

and, 



H{z) = {Z + Ij^Vn + ^(^ + l)AVn^ (44) 

respectively. Other physical magnitudes of observational interest are the scalar field equation 
of state (and barotropic index): 

and the deceleration parameter: q = — (1 + H/H'^); 

g = -1 + AV2 + (37/2)0^. (46) 

3.1.1 Relationship among frames 

Now it is useful to give the main physical magnitudes of observational interest also in the 
Jordan frame, as related with the EF magnitudes. Recall that dt — e 

a = e-^/^V^a = a" (see the definition for n in Eq. (42)), and f^b = e^^l^^T^^ ^ Pm ^ 

Other important relationships are the Hubble parameter: 

H = ne^/^'^^^H = nd^^^^^^W, (47) 
the dimensionless density parameter, 

= (l/n^)^^, (48) 
the dimensionless scalar field density parameter,^ 

= 1 - (1 V)0„, (49) 

the scalar field equation of state (and the barotropic index): 

*In the Jordan frame we define fli = Pi/ PcrUicai, where pcriUcai = 'ie'^H^, so that the Friedmann constrain 
equation (26) can be written as: 1 = Qm + ^4>- 
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. , W3)(A-l/v/e)-|(7/Ve)0.. .... 

^0 = -1 + 70, !<!> = A ^2 _ o ' ^^^^ 

and the deceleration parameter, 

q = (l/n)[g-+ (1 - n)(l - AV2) - (n/2)A/V^]. (51) 

This class of solutions depends (in principle) on 5 parameters: A, ^, po, Vq, and 7. The 
constants A and B are given in terms of the aforementioned parameters. 



3.2 The space of parameters 

In this subsection we will reduce the space of parameters by properly choosing the normaliza- 
tion for the Hubble parameter and by using some observational facts. In the Einstein frame 
the Priedmann constrain equation (29) means that Qm < 1- Besides, the solution Q,^ — 1 
is not stable so we are left with Clm < 1- The same is true in the Jordan frame, where 
the Fricdmann constrain equation (26) means that < Qm < 1- On the one hand, Qm is a 
maximum at 2; = 00 (the same is true for flm,)- Then, we require that ^2^(00) = l — e, where, 
as before, £ is a small number. This, in turn, yields: fim(oo) = n^(l — e) In consequence, 
Po,7 = 3n^A(l — e). On the other hand, according to the observations, at present (in the JF) 
^m(O) = (l/n^)r2m(0) = 1/3. This last equality implies that po,7 = n'^{A + B). Combining 
of the above constrains yields B/A = 2 — 3s. In consequence, Eq. (43) and Eq. (44) can be 
written in the following way: 

(z + 

Q^(z) =n^(l-6)- ^ ,,,, , (52) 

^ ^ ^ ^(^ + l)('5^-^')/" + 2-3£ ^ ^ 



and. 



\/P0,7 



n 



(^+l)^^/" + (2-3£)(^+l)^Vn 

^ , (53) 



\ 3(1 



respectively. The space of parameters (7, A, e, ^0,7) can be further reduced by choosing the 
normalization in which H{0) = 1 =^ -^(0) = ^/n =^ Po,7 = 1- Iii this case we are left with 
a three-parametric class of solutions spanned by the free parameters (7, A, e). In section 4 we 
will further constrain the free parameters A and e by means of comparison with supernovae 
observational data. 



3.3 Existence and Stability of the Solutions 

The stability of the system of dynamical equations (29-31) has been already studied in ([19]) 
for a single exponential of the kind V = Vqc^^^ . This study has been performed in the 
Einstein frame, however, the existence and stability of critical points already defined in the 
Einstein frame, remain unchanged under a conformal transformation back to the Jordan 
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frame [19]. The system (29-31) can be arranged in the form of a plane-autonomous system 
if one introduces the variables: 



In this case one obtains [19]: 

/3 3 



and, 



y' = -\l^gxy + ^7/(2x' + 7(1 - - y^)), (55) 



x' ^-3x + J ^gy' + ^x{2x'' + 7(1 - - y')) + W{1 - x' - y^). (56) 



Following our previous result in section 2.C, that the linear relationship between the time 
derivative of the scalar field and the Hubble parameter leads to solutions being attractors of 
two possible kinds: scalar field dominated and scaling solutions, then one expects that the 
relationship (34) leads to the same result in the present situation. 

The system described by (55) and (56) can have up to five fixed points depending on the 
values of 51,7 and W [19]. However, when combined with the relationship (34), that leads 
to the fixed value: x — A/\/6, it leaves us with two possibilities (see the potential in Eq. 



(35)): 1) Either g^^ = 1, and 7^ = AV3, or 2) c/2 = = 37/A + J2/3W =^ Q^, = 



A(A - y'2/3W^)/37, and % = 7A/(A - ^2/3W). 

In the first case, in the Einstein frame, we have a scalar field dominated solution, mean- 
while, in the second case we have a scaling solution. Solution 1) exists for A^ < 6 and is 
stable when W - \/W'^ + I87 < ^eA < FT + ^W'^ + I87. Meanwhile, solution 2) exists 
for A > 37iy/(^^^ - 1) and either A > ^fl or A < ^fl . It is stable 

whenever it exists [19]. 



4 Observational Testing 

Although we can not talk literally about observational testing of the solutions found (we 
point out that the observational facts are not yet conclusive in some cases and are not 
accurate enough in others), in this section we shall study whether our solutions agree with 
some observational constrains that are more or less well established. The main observational 

facts we consider are the following[2]: 

1. - At present [z = 0) the expansion is accelerated {q{0) < 0). 

2. - The accelerated expansion is a relatively recent phenomenon. Observations point to 
a decelerated phase of the cosmic evolution at redshift z — 1.7. There is agreement in that 
transition from decelerated into accelerated expansion occurred at a 2; 0.5 [24]. 
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3. - The equation of state for the scalar field at present a;</,(a;<^) ~ —1 (it behaves hke a 
cosmological constant). With a 95% confidence limit uj < —0.6 [25]. 

4. - Although, at present, both the scalar (quintessence) field and the ordinary matter have 
similar contributions in the energy content of the universe (0^(0) = 1/3 =^ f2<^(0) = 2/3), 
in the past, the ordinary matter dominated the cosmic evolution, ^ meanwhile, in the future, 
the quintessence field will dominate (it already dominates) and will, consequently, determine 
the destiny of the cosmic evolution. 

Now we proceed to "observationally" test the solutions found in the cases studied in the 
former sections. 

4.1 Minimally coupled quintessence 

In the case of normal quintessence model studied in section 2, our solution depends on two 
parameters e and k (we fix 7 = 1, meaning cold dark matter dominance at present). To 
constrain the values of these parameters we use the following facts: 

1. Corasaniti and Copeland [26] found that the determination of the third peak in the 
BOOMERANG data limits the value today of equation of state —0.93 < id^{0) < —1 

2. As stated before [27], nucleosynthesis predictions claims that at 95% confidence level. 
%{lMeV c^z = 10^°) < 0.45 

3. During galaxy formation epoch [28], around z ~ 2 — 4, the value of quintessence 
density parameter is Q,^ < 0.5 

Using these constrains we bound the parameter space (fc, e) by means of a computing 
code. These parameters range as < e < 0.045 and k > 0. 

Now we select a pair from the above interval in order to compare our model with Super- 
nova data. We first compare our model with the theoretical predictions of the Lambda Cold 
Dark Matter model {AC DM) for modulus distance. Subsequently, we follow the procedure 
used in Ref. [29] to further test our model. However, this procedure, which implies direct 
comparison of our model with observations of effective apparent magnitudes of SNla, can not 
be used to further constrain the space of parameters, since the free parameters in our model 
are not much sensible to these observations. For instance if we consider e — 0.01, the 
distribution has a minimum in mo direction^" at mo — 23.97 (x^ = 62.5996 or x^/dof — 1.18 
per degree of freedom) , however it has no minimum in the k direction, meaning that the 
parameter k can take any positive value. In fig. 1, is plotted as function of the free 
parameters k and mo (we chose e = 0.01, k could be any value in the physically meaningful 
range) . Although the model fulfils the observational requirements 1 — 4 given at the begin- 
ning of this section, we point out that other observations should be considered in order to 
further constrain the s])ace of ])arameters. 

sufficiently long matter dominated decelerated phase is needed for the observed structure to develop 
from the density inhomogeneities [15] 

^^mo is a parameter connected to the absolute magnitude and the Hubble parameter 
^^In the limit A; — > 00 the ACDM model is recovered. 
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4.2 Brans-Dicke quintessence 

In this case our solution depends on two parameters A and e (recall that we fix 7 to be 
unity). We apply exactly the same procedure that in the former subsection to constrain 

the parameter space. However, in this case we should differentiate among the two relevant 
frames: The JF and the EF, in which Brans-Dicke theory can be formulated. Wc treat both 
formulations separately, bearing in mind that none of them could be taken to be the one 
with real physical meaning. We recall that, given the uncertain character of the dark energy 
(it dominates the energy content of the universe) and dark matter (it dominates the matter 
density in the universe), it could be possible that the dark energy and the dark matter, on 
the one side, and the ordinary (barionic) matter on the other, are minimally coupled with 
respect to different metrics (see Ref.[19] and references therein). Therefore, both conformal 
formulations are of physical significance. 

When the aforementioned procedure is applied to the Einstein frame magnitudes of ob- 
servational interest, one finds that the physically meaningful region in parameter space is 
bounded by < 5 < 0.045 and < A < 0.37. Meanwhile, in the Jordan frame one obtains 
that the physically meaningful range is < £ < 0.045 and 0.03 < A < 0.38. 

It should be pointed out that, in all cases, the physically meaningful region in parameter 
space is chosen such that the main observational facts 1—4 explained at the beginning of 
the section are fulfilled. As an illustration, in fig. 2, we show the evolution of both dimen- 
sionless background and scalar field energy densities & z for the Einstein frame magnitudes 
VLra and VL^ respectively^^ (the values of the free parameters are taken from the physically 
meaningful region in parameter space). In fig. 3, the evolution of the equation of state in 
both formulations uj^ and o;^ & 2; is shown, meanwhile, in fig. 4, we plot the deceleration 
parameter in both frames q{z) and q{z) to show the transition redshift when the evolution 
turns from decelerated into accelerated {z* = 0.54 and z* = 0.56 in the EF and in the JF 
respectively). It is seen that the differences in these magnitudes are not significant as one 
goes from one frame into the other. In fig. 5, as an illustration, we plotted the effective 
apparent magnitude observed [30] and calculated according to our model with parameters 
{e = 0.01, A = 0.2) to show the good fitting (the results are within la expectation). 

We want to stress that, as in the former case (subsection 4.1), the free parameters of the 
model are not much sensible to the SNIa data so, the standard procedure used for instance 
in Ref. [29], is not suitable for further constraining of the space of parameters of the model. 
Other observational evidence could be considered for this purpose. 

5 Conclusions 

We have found a new parametric class of exact cosmological scaling solutions in gravity 
theory minimally coupled to a self-interacting scalar field and in Brans-Dicke theory with a 
self-interacting BD scalar, by assuming a linear relationship between the Hubble expansion 

^^The Jordan frame magnitudes ^lm and fi^ do not differ significantly from the corresponding EF ones. 
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parameter and the time derivative of the scalar field. This relationship is suggested by an 
implicit symmetry of the field equations. 

We point out that it is not necessary to make any a priori assumptions about the func- 
tional form of the self-interaction potential or about the scale factor behavior. These are 
obtained as outputs of the assumed relationship once one integrates the field equations ex- 
plicitly. As a result a peculiar class of potentials arise: a combination of exponentials with 
some dependence on the barotropic index of the ordinary matter fiuid. This kind of potential 
is of interest for models of dark energy. Recall that a single exponential potential can lead to 
one of the following scaling attractor solutions [12]: 1) a late time attractor where the scalar 
field mimics the evolution of the background barotropic fiuid, and 2) the late time attractor 
is the scalar field dominated solution. In consequence, a combination of exponentials should 
allow for a scenario in which the universe evolves through a radiation-matter regime (attrac- 
tor 1) and at some recent epoch evolves into the scalar field dominated regime (attractor 
2)- 

By inspection of the stability and existence of the solutions found, we show that the 
assumed linear relationship between the Hubble parameter and the time derivative of the 
scalar field can lead always (both in standard minimally coupled and in Brans-Dicke models 
of dark energy) to attractor scaling solutions. Otherwise, the assumed relationship acts as a 
selection principle: only two of five critical points in phase space agree with it. 

In all cases, the relevance of the solutions found is outlined and the space of parame- 
ters is constrained with the help of a computed code that considers some well established 
observational facts. The models explored here are tested with the SNla observational data. 
Agreement with observations is always within la expectation. It is pointed out, however, 
that the procedure of Ref. [29], can not be used to further constrain the parameter space, 
since the free parameters in the models of interest are not sensible to the luminosity distance 
SNIa observational data. Nevertheless, the scahng solutions found in both cases, meet the 
main observational facts considered in this paper with acceptable accuracy (points 1-4 at 
the beginning of section 4.). It is also pointed out that other observational data should be 
considered to further constrain the space of parameters. 

In the case we deal with minimally coupled quintessence, the model explored in this 
paper is a particular case of the one in Ref. [12]. However, in that paper the authors did not 
consider any exact solutions at all but their aim was to numerically explore the model. In the 
Brans-Dicke far as we know, the combination of exponentials has not been explored 

yet, although the stability and existence of scaling solutions in Brans-Dicke quintessence 
with a single exponential potential has been studied, for instance, in Ref. [19]. 

We thank our colleagues Rolando Cardenas, Diosdado Villegas and Carlos R. Fadragas 
for helpful comments and Indrajit Chakrabarty for pointing to us Ref. [19]. We acknowledge 
the MES of Cuba by financial support of this research. 
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Figure 1: The function is plotted as function of the free parameters k and mo (we chose 
£ = 0.01. As seen k could be any value in the physically meaningful range so, SNIa luminosity 
observations do not allow for further constrain of the parameter space. Other observations 
could be considered for this purpose. 



Figure 2: We show the evolution of both dimensionless background and scalar field energy 
densities & z for the Einstein frame magnitudes and Q^. The parameters A = 0.2 and 
e — 0.01 have been chosen so that the SNIa data is reproduced within la. Equality of matter 
and quintessence energy density occurs approximately at z 0.3 — 0.4. 



Figure 3: The evolution of the equation of state in both formulations ou^ and ou^ &i zis shown. 
The equation of state for the quintessence field in the Einstein frame is never positive. 



Figure 4: We plot the deceleration parameter in both frames q{z) and q{z) to show the 
transition redshift when the evolution turns from decelerated into accelerated: z* = 0.54 — 
0.56. It is seen that the differences in these magnitudes are not significant as one goes from 
one frame into the other. 



Figure 5: As an illustration, we plotted the effective apparent magnitude observed [30] and 
calculated according to our model with parameters (e — 0.01, A = 0.2) to show the good 
fitting. The results are within la expectation. 
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